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FOREWORD 


This final report, submitted to RCA Government Services in fulfillment of 
contract G788563-JA2, discusses the feasibility of a gas grain simulation 
facility for the Space Station. 

This work was conducted by the Planetary and Applied Sciences Unit in the 
Payloads and Sensors Section of Martin Marietta Space Systems. We are 
indebted to Ms. Judith C. Powelson for contributions to two sections of this 
report. We also wish to acknowledge many valuable discussions and helpful 
ideas by the Technical Monitor, Dr. Guy Fogleman. 
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INTRODUCTION 


This report discusses various physics aspects of gas/grain experiments 
in microgravity. Its primary purpose is to elucidate the problems that 
must be dealt with and in many cases to assign values to the factors 
involved. This work is intended as a preliminary survey-type approach 
and must be followed up in each individual specific experiment to 
confirm the feasibility of the design. The report is organized as 
follows: 

(a) forces on individual particles, 

(b) levitation and particle handling methods, 

(c) levitation induced coagulation forces, 

(d) feasibility of long duration experiments, 

(e) feasibility of classes of experiments, and 

(f) summary and recommendations. 
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2.0 


FORCES ON INDIVIDUAL PARTICLES 


This section covers forces that act on individual particles. A later 
section covers the net effect of levitation-induced forces that tend to 
coagulate the particles. 

2.1 INTERPARTICLE FORCES 

Interparticle forces are due to Van der Waals attraction, electrostatic 
attraction, surface tension, and mutual shielding from the gas. 

2.1.1 VAN DER WAALS FORCES 

Van der Waals forces are defined as those that form the basis for the 
constant "a" in Van der Waals equation of state for a gas. They 
include forces between molecules having dipoles and quadrupoles, forces 
resulting from the polarization of molecules by static dipole or 
quadrupole fields of other molecules, and the nonpolar Van der Waals 
forces. Except when the dipole moment is large and the polarizability 
is small, the nonpolar forces (called dispersion forces) will be the 
largest (Corn, 1966). 

The nonpolar Van der Waals forces between two particles depend on the 
electrodynamic properties of their respective media. Since there is a 
finite time of propagation of electromagnetic waves between the two 
particles, the effects are retarded and the potential is reduced (a 
very readable explanation of retarded potentials is given by Spruch, 
1986). Modern Van der Waals theory was developed by Lifshitz (Marlow, 
1980). For separation distance very much less than the wavelength of 
the interaction, Lifshitz theory reduces approximately to the 
(non-retarded) London Van der Waals dispersion forces as applied to 
condensed matter by Hamaker. Therefore for two spheres the attractive 
force for distance less than 1 nm is (Zimon, 1980) 


A 3i 02 

6 (a, + 82) ' 


( 2 . 1 . 1 - 1 ) 


where 

A 

aj , &2 
d 


Hamaker constant for the interacting media (J), 
= radii of spheres (m), and 
distance between surfaces (m). 


For separation distances very much greater than the wavelength of the 
interaction (>50 nm), Lifshitz theory reduces to the (fully retarded) 
theory (Marlow, 1980). Therefore for large d (Zimon, 1980) 


4 B a, 3; 

3 (a, + a^) ' 


( 2 . 1 . 1 - 2 ) 


where B (J-m) is an interaction constant. 

In aerosol microphysics, Van der Waals forces are important for 
particle separations less than a radius (Marlow, 1980). For d << a 
(radius of smaller sphere), an approximate expression for the partially 
retarded interaction potential energy between two smooth spheres is 
(Gregory, 1981) 


V = 


A a, 3; 

6 d (a, + 82 ) 



(2.1. 1-3) 


where 

X = characteristic wavelength of the electromagnetic interaction 
b ~ 5.32 . 

In Eq. 2. 1.1-3 the term in brackets accounts for retardation. The 
Hamaker constant may be computed from dielectric data for the two 
particles and the gas (Hough and White, 1980). Some theoretically 
computed values of the Hamaker constant for particles in air are listed 
in Table 2. 1.1-1. Some experimental values of the Hamaker constant for 
particles in a vacuum are 4.7 x 10“^^ J for borate and 2.2 x 
10“^^ J for quartz (Corn, 1966). 

For spherical particles with rough surfaces, the right hand side of Eq. 
2. 1.1-3 should be multiplied by the factor d/s, where d is now defined 
to be "the distance between outermost peaks of surface irregularities 
and s is the separation between the mean profile lines" (Czarnecki, 
1984) (See Fig. 2. 1.1-1). 
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Figure 2. 1.1-1 R ough Spherical Particles 


For nonspherical particles such as flat plates, rectangular rods, and 
cylinders, particle orientation has an important effect on 
interparticle forces. These particles usually orient themselves such 
that the largest plane areas are opposed, which minimizes the potential 
energy (Hidy, 1984), 

Even in a vacuum, the calculation of Van der Waals forces is at best 
approximate because of uncertainties in composition, size, shape, 
surface roughness, orientation, and spacing. In a gas, the dispersion 
force may be overwhelmed by adhesion forces due to electrostatic charge 
or absorbed molecules, in which case the dispersion force calculation 
is valuable only as an order of magnitude estimate (Corn, 1966). 

Van der Waals forces scale with particle radius while gravitational 
forces scale with radius cubed. Thus Van der Wals forces become less 
important with increasing particle size (Corn, 1966). 

Temperature affects Van der Waals forces. However, temperature effects 
that are confined to the oscillator modes of a particle may be 
neglected for (hc/kT) » diameter where T is the particle temperature 
(Marlow, 1980). 

Table 2. 1.1-1 Theoretical Values of the Hamaker Constant for Particles in 

Air (from Hough and White, 1980). 


Material of Particle 

Crystalline quartz 
Fused quartz 
Fused silica 
Calcite 

Calcium fluoride 
Sapphire 

Poly (methylmethacrylate) 
Poly (vinylchloride) 
Polystyrene 
Polystyrene 
Poly (isoprene) 

Poly (tetrafluoroethylene) 
"Teflon FEP" ^b) 


Hamaker Constant (xlO~^®J) 


8.83 

6.50 

6.55 

10.1 

7.2 

15.6 

7.11 

7.78 

6,58 

6.37 

5.99 

3.80 

2.75 


(a) variation in refractive index between sources 

(b) using infrared data for poly (tetrafluoroethylene) 
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2.1.2 ELECTROSTATIC INTERPARTICLE FORCES 

Particles in a cloud will be charged. Even if the cloud is not charged 
initially it will become charged by diffusion of ions present in the 
gas and by contact charging. Contact charging is not due to friction 
but results from the making and breaking of contacts. The cloud will 
contain both negatively and positively charged particles with the 
equilibrium charge distribution being approximately Gaussian. The 
charges increase slightly less rapidly than the square of the particle 
diameter (Corn, 1966). 

The initial charge on liquid droplets may be due to spray 
electrification. At an interface between a gas and certain liquid 
substances, electric dipoles are created with positive polarity inward 
and negative outward. When the surface is disrupted to form fine 
droplets, charged particles result (Corn, 1966). 

The origin of initial charge on solid particles is not well 
understood. It may be complicated by the effect of surface impurities, 
films of moisture, formation of Helmholtz double layers or by 
electrolytic effects. Electrification can also be affected by the 
state of the surface and heat and mechanical stresses of the material. 
Also, certain materials may transfer ions across surfaces after initial 
contact (Corn, 1966). 

Contact potential differences may also arise due to the surface effects 
of dissimilar materials. Their electrostatic binding force is 
proportional to the particle dimension to the two-thirds power 
(Barengoltz and Edgars, 1975). 

There is much experimental evidence for the great importance of 
electrostatic forces in the adhesion of particles. On the other hand 
crude calculations show that for the electrostatic forces to exceed the 
Van der Waals forces the charges would have to be much larger than 
observed and would have to be greatly concentrated near points of 
contact (Corn, 1966). 

A more recent book (Hidy, 1984) does not quantitatively improve much on 
the above references to Corn. Hidy concludes his section on this 
subject with a call for more study:” . . . future work on the 
fundamental nature of adhesion must incorporate studies of electrical 
interactions between solids, as well as intermolecular interactions.” 
Hidy's only theoretical reference on this subject is Derjaguin and 
Smilga (1967), no references to which have appeared in the Science 
Citation Index since 1982. 

A computerized search of the last 10 years of the IEEE (including 
Physics) Abstracts and the Aerospace Abstracts for the keywords 
["electrostatic?” with "interparticle” with ("force?” or "attraction" 
or "coagulation”)] yielded 20 references none of which is relevant to 
the above discussion. 
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2.1.3 SURFACE TENSION 


For a wet or dirty surface, the capillary effects of surface films may 
produce an additional attractive force. Capillary condensation causes 
a film of water to form between the particles. For a relative humidity 
of lOOX and for two equal sized spheres in contact, the adhesive force 
due to surface tension is 


F « 2n Rp o , 


(2.1. 3-1) 


where Rp is the radius of the particles and a is the surface 
tension. Eq. 2. 1.3-1 has been verified by direct experimental 
measurement. The force will be less for a lower relative humidity. 

The adhesion due to surface films may exceed that resulting from either 
intermolecular (Van der Waals) or electrical forces (Hidy, 1984). 

2.1.4 MUTUAL SHIELDING FROM THE GAS 
The Knudsen number is 


Kn = (2. 1.4-1) 


where 

Xg = mean free path of the gas (m), and 
Rp= particle radius (m). 

Table 2. 1.4-1 lists the mean free path of nitrogen at 20*C. 

As Kn — 0 continuum theory applies, whereas as Kn ►-the dynamics 

is characterized by individual gas molecules colliding with the 
particle (Hidy, 1984) and is called the free-molecular regime (Marlow, 
1980). For 1 s Kn $ 10 the particle is in the transition regime. For 
two nearby particles in the transition regime, gas-particle collisions 
are less frequent between the particles and the particles experience an 
effective attractive force (Marlow, 1980). 
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Table 2. 1.4-1 Mean Free Path of Nitrogen at 20®C 
(from CRC, 1982) 


Pressure 
(mm Hg) 

1 

0.1 

0.01 

0.001 


Mean Free Path 
(m) 

4.5x10“^ 

4.5x10-'^ 

4.5x10-3 

4.5x10-2 


Note: 1 mm Hg = 1 torr = 1.33 millibar. 


I 
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2.1.5 OTHER FACTORS 


An external electric field has an immediate effect for example by 
producing aggregate threads of particles. It may also have a lasting 
effect in that polarization of particles has been found to persist 
after the field was shut off. Other factors influencing adhesion 
include (a) chemical and electrical properties of the materials, (b) 
surface contamination, nature of the particle surface contact (particle 
size, shape, roughness and state of aggregation), (c) temperature, and 
(d) time of contact (Corn, 1966). 

2.2 DRAG FORCE 

Particle motion through a gas may be characterized by the Reynolds 
number Re and the Knudsen number Kn. The Reynolds number is 


Re = 



( 2 . 2 - 1 ) 


where 

Rp = particle radius (m), 

Vp = particle velocity relative to gas (m/s), 
Pg = mass density of gas (kg/m^), and 
y® = viscosity of gas (kg/m-s). 


The drag force acting on a particle is 



(2.2-2) 


9 


where B is the mobility. For Kn— *-0, Re « 1 , the mobility of a 
spherical particle is 


6 n |J, Rp ' 

For K„ — 0, Re < 1 , 

F = ^ p, n Rp^ C Vp* , and 
C = (12/Re)(1 + 0.250 Re"'’). 


For K„ 0, Re < 1 , 


B = 
A = 


6 n n, Rp 


, and 


1 + 1.257 Kn + 0.400 Kn exp (-1.10 Kn ’). 


(2.2-3) 

(2.2-4) 

(2.2-5) 

( 2 . 2 - 6 ) 

(2.2-7) 


Nonspherical particles are included by use of dynamic shape factors 
(Hidy, 1984). 

2.3 BROWNIAN MOTION 

The mean square displacement of a particle due to Brownian motion is 
(Spumy, 1986) 


Ax^ = 2tkTB, 


(2.3-1) 


where 

t = time(s), 

k = Boltzmann constant (1.3804 x 10~23 j/k), 
T = gas temperature (K), and 
B = mobility of the particle. 
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The number of particles striking the entire inner surface per unit time 
of a sphere of radius a is (Hidy, 1984) 


8a nDp N. 



-n^ D 



(2.3-2) 


where 

N, = initial particle concentration far from the wall, and 
Dp = BkT. 

It is useful to compare Brownian motion with that due to the terminal 
velocity at 10~5 g. The particle is placed at the center of a 
spherical chamber of radius x. At terminal velocity the drag force of 
Eq. 2.2-2 equals m(l0“5 g). Assuming that the terminal velocity is 
reached instantaneously, the time until the particle hits the spherical 
wall is 


9 tip X 

2Rp^ Pp A (10 »g) ' 


(2.3-3) 


where is the mass density of the particle and A is given by Eq. 
2.2-7. The time for the particle to reach the wall by Brownian motion 
is 


U 


3 n |i^ R„ 
kT A 


(2.3-4) 


For a chamber of radius 0.5 m, for nitrogen gas at 300 K, and for 
particles of mass density 1 gm/cm^, particle radii corresponding to 
various values of tg/tg are listed in Table 2.3-1. For a chamber 
of radius 0.5 m, nitrogen gas at 300 K, and a particle of mass density 
1 gm/cm^, these times are evaluated for various particle radii and 
gas pressures in Tables 2.3-2. and 2.3^3. From Table 2.3-2 it may be 
noted that at a pressure of 1 mm Hg, Brownian motion causes a particle 
of radius 1 nm to reach the wall in only two minutes. In particular, 
the nanometer sized grains that are nucleated at a pressure of 0.2 mb 
(0.15 mm Hg) in the Titan aerosol experiment (McKay et al., 1986) would 
reach the wall by Brownian motion in 20 seconds [if they were not 
formed in some kind of levitation well and if they did not first 
coagulate (coagulation depends on particle number density, which was 
not specified)]. 
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Table 2.3-1 Radii for Various Ratios of Times, until Particle Hits Wall 
by Brownian Motion and 10”^g. 


‘B/t| 

Rp(ym) 

0.1 

0.16 

1 

0.34 

10 

0.74 


Conditions: 

T = 300 K 
= 1 g/cm^ 
X = 0.5 m 
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Table 2.3-2 Time Until Particle Hits Wall 

Due to 10”5g and Brownian Motion 
(Gas Pressure = 760 mm Hg) 


time (s) 


Rp(ra) 

Kn 

A 

10-5 g 

Brownian 

10-y 

59 

98 

4.2xl0l^ 

1 .0 xl 05 

10-8 

5.9 

10 

4.1xl0ll 

1.0 xl 07 

10-7 

0.59 

1.8 

2.3x10^0 

5.6 x 108 

10-8 

0.059 

1.1 

3.7 x 108 

9.2 x 109 

10-5 

0.0059 

1.0 

4.1 x 108 

1 .OxlOll 

10-^ 

5.9x10-"^ 

1,0 

4.1x10^ 

1. 0x10^2 

10-3 

5.9 x 10-5 

1.0 

4.1 x 102 

1.0xl0l3 


Conditions: 

Chamber radius : 
Gas: 

Temperature: 
Particle density: 


0.5 m 
Nitrogen 
300 K 
1 g/cm^ 
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Table 2.3-3 Time Until Particle Hits Wall 

Due to 10“5g and Brownian Motion 
(Gas Pressure = 1 mm Hg) 


time (s) 


Rp(m) 

Kn 

A 

10”5g 

Brownian 

10-y 

4.5x10'^ 

7.5x10^ 

5.4x10^ 

1.4x10^ 

10-8 

4.5x108 

7.5x108 

5.4x10® 

1.4x10'^ 

10-'^ 

4.5x102 

7.5x102 

5.4x10^ 

1.4x10® 

10-8 

45 

75 

5.4x10® 

1.4x10® 

10-5 

4.5 

8.1 

5.0x105 

1.2xl0l0 

10-^ 

0.45 

1.6 

2.6x10^ 

6.3xl0ll 

10-8 

0.045 

1.1 

3.7x102 

9.2x10^2 


Conditions: 

Chamber radius: 
Gas: 

Temperature: 
Particle density: 


0.5 m 
Nitrogen 
300 K 
1 g/cm^ 
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2.4 EXTERNAL QUASISTATIC FIELDS 

The force due to external fields is 

r = rrip? + Qp (E* + ^ X §1, 


where 

nip = mass of the particle (kg), 
g = acceleration of gravity (m/s^), 
qp = charge on the particle (C), 

E = electric field (V/m), 

Vp = particle velocity (m/s), and 
B = magnetic field (T). 
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ALTERNATING MAGNETIC FIELD 


An alternating magnetic field will generate eddy currents in a 
spherical particle. The eddy currents will be strongest at the 
surface. By Lenz's law, the currents will act to oppose the field and 
tend to exclude it from the inside of the particle. For a field Bq 
that varies slowly in the axial direction, the axial force is 


F 


- n 


H. 


V B,^ 


( 2 . 5 - 1 ) 


where a is the radius of the particle (Frost and Change, 1982). 

2.6 G- JITTER 

Random accelerations that affect a particle are called g-jitter and are 
due to crew motion, docking, equipment operation, etc. They may be 
approximated by a time varying acceleration (McKay et al., 1986) 


a, = So sin (tot), (2.6-1) 


where 

aQ = amplitude of the acceleration (m/s^), and 
oj = angular frequency (s“1). 

For Kn •+■ 00 and the maximum displacement Xq » it is called 
ballistic g-jitter and the particle behaves as if it were in a vacuum. 
For Xq drag forces must be considered: for ECn o» it is 

called Knudsen g-jitter and for Kn -► 0 it is called Stokes g-jitter 
(McKay et al., 1986). 

2.7 PHORESIS 

If an aerosol particle is in a temperature gradient, light, or 
vapor-concentration gradient, the Brownian motion is biased toward one 
direction and this is called a phoresis. Thermophoresis, 
photophoresis, and diffusiophoresis will occur, respectively (Spumy, 
1986). 
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2.7.1 THERMOPHORESIS 


The force due to a temperature gradient on a spherical particle is 
(Hidy, 1984) 

f ^ Rp" & VTexp (^\ . 12.7.1-1) 

15 V, \\J 

for 0.25 < A,/Rp < oo, Ma < < 1, 
where 

v_ = mean thermal speed of the gas (m/s), 

o 

Xg = mean free path of the gas (m), 

= Mach number, 

M. = Vp/v„ (2.7. 1-2) 

T = 0.9 + 0.12 a„, + 0.21 cx,„ ^1 - (2.7. 1-3) 

(for monatomic gases), 

Ojj s accommodation coefficient for momentum, 

= accommodation coefficient for heating, 

kg = thermal conductivity of the gas (J/m-s-K), and 

kp = thermal conductivity of the particle (J/m-s-K). 

The accommodation coefficients vary from zero to unity and must be 
measured for each material combination (Hidy, 1984). 
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2.7.2 PHOTOPHORESIS 


Photophoresis is due to heating of the particle on one side by light. 
The gas on that side is heated and the particle moves in the opposite 
direction. In positive photophoresis, the particle moves in the same 
direction as the light, while in negative photophoresis it moves toward 
the light source (Preining, 1966). 

In a beam of light, a sphere may form a lens, causing the light to be 
refracted within the particle. There are two cases. For a weakly 
absorbing particle, the image of the light source is close to the dark 
side of the particle. A small weakly absorbing particle will be heated 
more on the dark side and will exhibit negative photophoresis, while a 
large particle (due to absorption) will be heated more on the lighted 
side and will exhibit positive photophoresis. A strongly absorbing 
particle will be heated more on the lighted side and will always 
exhibit positive photophoresis (Preining, 1966). 

Negative photophoresis may be observed by shining laser light on smoke 
particles (which are not totally opaque) in air (Ashkin, 1972). 

2.7.3 DIFFUSIOPHORESIS 

The diffusion force depends on the concentrations, gradients, molecular 
masses, and diffusion coefficients of the gas species involved and also 
on the cross section of the particle (Hidy, 1984). 

2.8 RADIATION PRESSURE 

Light carries momentum and thus exerts force on a particle when it is 
reflected, refracted, or absorbed by it. The momentum carried by a 
photon is 


P 



(2.8-1) 


where 

E = energy of the photon (J), and 
c = speed of light (2.998 x 10^ m/s). 
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An absorbed photon imparts all of its momentum to the particle- A 
photon elastically scattered by 180® imparts twice its momentum to the 
particle. A photon refracted or reflected at some other angle imparts 
momentum according to the laws of conservation of energy and momentum. 
For a light intensity gradient perpendicular to the direction of the 
light, a component of force is applied perpendicular to the direction 
of light. For highly reflective particles, the force is toward the 
region of lesser intensity. For transparent particles with index of 
refraction greater than that of the surrounding media, the force is 
toward the region of greater intensity (Ashkin, 1972). 

2.9 ACOUSTIC FORCE 

A standing sound wave produces through a nonlinear effect an average 
ambient pressure change of (Wang, 1986) 
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where 

^ = time average of square of instantaneous pressure (Pa^), 

p = mass density of gas (kg/m^), 

c = speed of sound in gas (m/s), and 

^ = time average of square of instantaneous gas velocity (m^/s^). 

When placed in a standing wave field, a sphere experiences a radiation 
pressure force of (Wang, 1986) 


\ 6 / \pc7 


F 


( 2 . 9 - 2 ) 


where 


Pl= pressure amplitude of the fundamental frequency (Pa), 

k = wave number (m~^), 

a = sphere radius (m), and 

X = position of center of sphere (m). 

Eq. 2.9-2 is the basis of acoustic levitation and has the following 
limitations: 

(a) the mass density of the sphere must exceed that of the gas; 

(b) the diameter of the sphere must be much less than the 
wave length; and 

(c) the diameter of the sphere must be greater than the mean free 
path of the gas (McKay et al., 1986), 

It is obvious from Table 2. 1.4-1 that acoustic levitation can not 
confine a nanometer sized particle that first nucleates from a few 
molecules of the product of a gaseous chemical reaction. Also, 

Appendix C shows that for frequencies very much less than MHz, acoustic 
levitation prevents Brownian motion of the levitated particle. 

It should be noted that the force given by Eg. 2.9-1 is proportional to 
the volume. The sphere also experiences the viscous drag force of 
acoustic streaming, which depends on the surface area of the sphere and 
is a complicated function of (a) the geometry of the transducer, (b) 
the schemes for the production of standing waves, and (c) the acoustic 
boundaries (Lee and Feng, 1982). The levitation becomes unstable when 
the streaming force becomes as high as 25% - 75% of the volumetric 
force. As the size of the particle is decreased, the ratio of the 
streaming force to the volumetric force increases. However the ratio 
may be reduced by using a higher frequency (Lee and Feng, 1982) (see 
Fig. 2.9-1). M. Barmatz of JPL (private communication via G. Fogleman 
of RCA, July, 1987) has estimated that streaming limits the minimum 
particle radius to 1% of the wavelength. 
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3.0 LEVITATION AND PARTICLE HANDLING METHODS 
The levitation methods considered are 

(a) none (chemically inactive walls), 

(b) light pressure, 

(c) radiometric, 

(d) acoustic, 

(e) electrostatic, 

(f) aerodynamic, and 

(g) electromagnetic. 

3.1 NO LEVITATION 

Many cloud type experiments do not require levitation but should use 
inert walls. This would avoid the artificial coagulation caused by all 
of the levitation methods. Many walled cloud experiments that are 
currently performed on Earth would be much better in space because 
particles could grow to much greater size before they precipitated due 
to the residual gravity. 

3.2 LEVITATION BY LIGHT PRESSURE 

Levitation by light pressure is possible for highly transparent 
particles and for highly reflective particles (Ashkin, 1972). 

Partially transparent (absorptive) particles would be heated to the 
point where melting or chemical changes would have to be considered and 
in a gas would be subject primarily to radiometric (photophoretic) 
forces. Highly reflective particles [reflectivity >> 98% (Ashkin, 
1972)] would not correspond to the cases of interest. Therefore the 
following will be restricted to highly transparent particles. 

In a microgravity environment two opposing beams would be used (see 
Fig. 3.2-1). The index of refraction of the particle is greater than 
that of the gas. The laser beam has a Gaussian intensity distribution, 
which causes the centering of a particle in the beam. The particle is 
also pushed with the beam, and the diffusive spreading of the 
oppositely directed beams produces a trap for the particle (Ashkin, 
1972). 

For sufficient gas pressure, the fluid damping (Ashkin, 1972) would be 
sufficient for stable confinement of many particles simultaneously. In 
a vacuum, active feedback of laser power might have to be used and thus 
confinement might be restricted to one particle. 




Figure 3.2-1 

Optical Bottle Refraction for Centering Small Transparent Particles (Ashkin, 1972) 


It is possible to confine nonspherical particles. However, they become 
oriented in the beam (Ashkin and Dziedzic, 1980) and this would prevent 
a realistic study of coagulation processes. 

It is also possible to achieve stable confinement by an alternating 
light trap (Ashkin and Dziedzic, 1985). However, this would not have 
any advantages and would require many more optical components. Another 
method is to use the strong light gradient forces at a single optical 
focus (Ashkin et al., 1986). However, this restricts the trap to a 
very small size and would thereby preclude multiple particle 
experiments. 

A small number of particles (2-20) might be confined by directed pulses 
of laser light. 

3.3 RADIOMETRIC LEVITATION 

For highly absorbing particles in a gas, the radiometric 
(photophoretic) forces greatly exceed those due to light pressure 
(Lewittes and Arnold, 1982). A stable trap has been demonstrated in 
which a TEMOl* (doughnut) mode laser shines upward and suspends highly 
absorbing particles against the force of gravity. The TEMOl* mode has 
an intensity minimum at the center and thus highly absorbing particles 
are confined both radially and vertically (see Fig. 3.3-1) (Lewittes 
and Arnold, 1982). It seems reasonable that two opposing TEMOl* lasers 
could create a stable trap in microgravity (see Fig. 3.3-2). This 
would certainly be true if feedback control of laser power were used to 
confine a single particle. However, even without feedback a cloud of 
particles could probably be confined. 

The situation with less highly absorbing particles is more complicated 
and depends on their size and the gas pressure. Small particles will 
be subject to negative photophoresis and will move toward the light 
source [but also perpendicular to the light in the direction of an 
intensity minimum (Lewittes and Arnold, 1982)]. A single small 
particle could be confined by applying feedback to the configuration of 
Fig. 3.3-2 (a laser's intensity would be increased to attract the 
particle toward it). There does not appear to be a way of confining a 
cloud of particles small enough to be subject to negative 
photophoresis. Also, for weak enough absorption and a very high or low 
gas pressure, the light pressure forces may predominate over the 
radiometric forces. This is because the radiometric force is maximized 
at the pressure at which the mean free path of the gas is equal to the 
diameter of the particle (Ashkin and Dziedzic, 1976). 
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Figure 3.3-1 

Radiometric Levitation of Highly Absorbing Particle at 1 g 
(Lewittes and Arnold, 1982) 
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Larger and less highly absorbing particles are subject to positive 
photophoresis and would behave similarly to the highly absorbing 
particles described above. This would be true if all particles were 
greater than a certain critical size. However, it does not seem 
possible to confine a cloud of small particles (subject to negative 
photophoresis) and then continue to confine them as some grow to sizes 
subject to positive photophoresis. This is particularly true since as 
the size of a given particle reaches that at which positive and 
negative photophoresis cancel the predominate radial force would be 
that of light pressure, which would force the particle to the outside 
of the TEMOl* laser beam. This effect may preclude radiometric 
confinement during nucleation and growth of even a single particle, 
although perhaps a careful procedure may circumvent this for a single 
particle. 

3.4 ACOUSTIC LEVITATION 

The acoustic approach is one of the most highly developed types for 
levitation in microgravity and on Earth. Here, the emphasis will be on 
the leading types of acoustic levitation for microgravity: single-axis 

interference and tri-axis. However, other types will be mentioned 
because they may fill special needs. 

3.4.1 SINGLE AXIS ACOUSTIC INTERFERENCE LEVITATION 

In single axis interference levitation, a small reflector 
(diameter - X .= sound wavelength) is placed an arbitrary distance 
from a fixed frequency ultrasonic source. The reflected sound wave 
interferes with the incident sound wave to produce a region of minimum 
acoustic potential energy at a distance of slightly more than A/4 from 
the reflector (see Fig. 3.4. 1-1) (Day and Ray, 1986). This potential 
well, which is deep in the axial direction and shallow in the radial 
direction, traps the particle (Day and Ray, 1986). The reflector and 
driver are made slightly convex to focus the sound (Naumann and 
Elleman, 1986). 

The walls of the chamber are designed to absorb sound (Naumann and 
Elleman, 1986) and plane surfaces are avoided (Whymark et al . , 1979). 
Careful choice of chamber dimensions avoids resonances (Naumann and 
Elleman, 1986). The gas absorption is selected carefully: if it is 

too high, the energy well will be too weak; if it is too low, 
reflections may be a problem (Whymark et al., 1979). 
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Figure 3. 4.1-1 

Single-Axis Interference Acoustic Levitator (Day and Ray, 1986) 
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There are a number of advantages with this technique. First, no tuned 
cavity is involved. Therefore as the gas density, composition, 
temperature, or the sample size changes, the equilibrium position 
simply moves slightly relative to the reflector (Naumann and Elleman, 
1986). High specimen cooling rates (30°C/s) are possible (Day and Ray, 
1986). Levitation has been demonstrated up to 1600*C (Whymark et al . , 
1979), although the levitation force diminishes approximately as the 
inverse 5/2 power of temperature (Naumann and Elleman, 1986). 

Another advantage is that the sample may be moved by moving the 
reflector. In fact, multiple samples may be manipulated by the use of 
multiple reflectors (see Fig. 3.4. 1-2 (a)). Liquid drops may be merged 
by tilting the reflectors until the potential wells overlap (see Fig. 
3.4. 1-2 (b)) (Whymark et al . , 1979). Access is good from all sides 
except that of the sound source. 

The primary disadvantage is the shallowness of the potential well in 
the radial direction. The radial restoring force is provided by the 
weak effect of the gas streaming past the sample (Bernoulli effect) 
(Naumann and Elleman, 1986). A sample that is well confined axially 
might be lost due to the transverse acceleration caused by astronaut 
motion. An associated disadvantage is that a liquid sample will not 
have a precisely spherical shape because the equipotential surfaces of 
the well are not spherical. In fact, large droplets of liquid having 
low surface tension (such as water) are broken up (Naumann and Elleman, 
1986). Also, there is no independent control over sample rotation. 

3.4.2 TRI-AXIS ACOUSTIC LEVITATION 

The tri-axis levitator is comprised of three perpendicular sound 
sources and a cavity which is tuned for the fundamental mode for each 
of three perpendicular directions (see Fig. 3. 4. 2-1) (Barmatz, 1982). 
The sample is positioned at the middle of the cavity (Barmatz, 1982). 

By using a rectangular chamber with a square cross section and by 
maintaining a 90® phase difference between the driver located on the 
equal sides, the sample can be rotated. Vibrations can be induced in a 
liquid drop by causing a sinusoidal amplitude modulation in one of the 
drivers (Wang, 1983). 

The primary advantage of the tri-axis levitator is that the acoustic 
well has spherical equipotential surfaces thus trapping the sample 
equally strongly from all directions and causing a liquid sample to 
remain spherical (Naumann and Elleman, 1986). 

A disadvantage is that a change in gas temperature (Naumann and 
Elleman, 1986), composition or density requires adjustment of the 
frequency. Although changes in sample size also require frequency 
adjustment (Barmatz et al . , 1983), the sample sizes of interest in 
gas-grain simulations will in most cases be too small to have any 
effect . 
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Figure 3.4. 1-2 

Sample Manipulation Using Multiple Reflectors 
(Whymark, etal., 1986) 


30 




Figure 3. 4.2-1 Triple-Axis Acoustic Levitator (Barmatz, 1982) 
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Another disadvantage is that only one sample can be levitated at a 
time. While the cavity walls may be constructed of plexiglass to 
enable visual inspection (Wang, 198A), the walls restrict mechanical 
access. Also, the gas temperature of the tri-axis system is currently 
limited to 600°C (Day and Ray, 1986). A limitation on the rate of 
temperature decrease (at constant position) to a few degrees per minute 
(Naumann and Elleman, 1986) would probably not affect most gas-grain 
experiments. 

3.4.3 OTHER CONFIGURATIONS FOR ACOUSTIC LEVITATION 
Some other configurations are the following: 

(a) hemispherical focusing radiator, 

(b) siren, 

(c) tetrahedral arrangement of sound sources, 

(d) ring-type radiator, and 

(e) single-mode levitator. 

Items (a), (b) and (e) were designed for use in 1 g, but might be 
useful for microgravity. 

The hemispherical focusing radiator uses a spherical concave reflector 
placed one half of a wavelength beyond the focus of the radiator. The 
back surface of the radiator is covered with 130 piezo-electric 
transducers. The advantage of this configuration is that it 
concentrates a large sound power in a small area. Samples of 
submillimeter size have been levitated. One disadvantage is that it 
causes the sample to rotate (Lee and Feng, 1982). 

The siren was designed to levitate large objects on Earth (Gatnmel et 
al., 1982). Its high power is not necessary for the levitation of 
small grains in microgravity and may cause noise pollution on the Space 
Station. 

The tetrahedral arrangement of sources does not require any reflectors, 
but still offers somewhat limited access (see Fig. 3. 4. 3-1) (Hatano et 
al., 1982). 

The focused ring source was designed to produce easy access around the 
object by using a single radiator (see Fig. 3. 4. 3-2). A piezoelectric 
- ceramic transducer drives the ring-type radiator through a solid 
horn. The ring is uniformly vibrated in the radial direction. Since 
the inside of the ring is tapered at an angle of 20®, the sound is 
focused a short distance from the ring (Hatano et al., 1982). 

The single-mode acoustic levitator uses a single non-planar mode 
generated from a single transducer. It uses reflections from all of 
the walls and an adjustable plunger controls the chamber length (see 
Fig. 3.4. 3-3). The "inherent phase coherence of the orthogonal force 
components automatically assures ..." translational and rotational 
stability. By exciting more than one mode the sample can be rotated in 
a controlled manner (Barmatz, 1987). 
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Figure 3. 4.3-1 

Tetrahedral Arrangement of Sound Sources 
(Hatano, et al., 1982) 
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Figure 3A.3-3 Single • Mode Levitator (Barmatz, 1987) 
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3.5 


ELECTROSTATIC LEVITATION 


Ernshaw’s theorem states that stable confinement with static electric 
fields is impossible (Rhim et al., 1985). Therefore, it is necessary 
to use either a feedback system or alternating electric fields. 

Electrostatic levitation with feedback will work for particles of any 
size and has been demonstrated in microgravity. Levitation with 
alternating electric fields has been demonstrated with small particles 
in 1 g and may work in microgravity. 

Electrostatic levitation works in a gas or vacuum with liquid or solid 
particles that are conducting or nonconducting. 

A dielectric particle having no net charge may be levitated in a 
nonuniform electric field. Appendix B shows that this is probably 
infeasible for gas-grain applications. 

3.5.1 ELECTROSTATIC LEVITATION WITH FEEDBACK 

The best electrostatic levitation method for large objects appears to 
be the tetrahedral positioner, in which spherical electrodes are 
located at the verticles of a tetrahedron (see Fig. 3. 5.1-1). Two 
cameras sense the position of the object and a proportional- 
integral-differential feedback algorithm is used to center it. The 
advantage of the tetrahedral positioner over the single-axis (eg. 
ring-type with feedback) positioner is that the damping (derivative 
term) is applied to oscillations in all directions. Since the 
development of the tetrahedral positioner emphasized application to 
materials processing, any problems with application to very small 
particles were not detailed. (Rhim et al., 1985). 

Initial charging of a (solid) object to be levitated may be by contact 
with one of the electrodes. For levitation for long periods and 
especially at high temperatures, additional charging may be provided by 
a pulse-charging system (Makin et al., 1986). 

3.5.2 LEVITATION WITH ALTERNATING ELECTRIC FIELDS 

A very convenient method for levitating and studying an aerosol 
particle in 1 g is through the combination of static and alternating 
electric fields shown in Fig. 3. 5. 2-1. The static voltage 
balances gravity. For a sufficiently high frequency, alternating 
voltage applied to the central ring the particle will remain 
stationary. No feedback system is required (Straubel and Straubel, 
1986). The reference did not make clear if the particle would remain 
stationary in microgravity with 0. 

3.6 AERODYNAMIC LEVITATION 

The following types of aerodynamic levitation are orientation- 
independent and therefore may be used in microgravity: 
constricted-tube gas-flow, triaxial sonic pump, and triaxial molecular 
beam. Plasma levitation may also be configured for microgravity. 
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Figure 3. 5, 1-1 Tetrahedral Electrostatic Positioner 


O 




37 



I I Pl«te 


Figure 3. 5. 2-1 Levitation with Static and Alternating Electric Fields 



3.6.1 CONSTRICTED-TUBE GAS-FLOW LEVITATION 


The tube contains identical converging and diverging regions on the 
upstream and downstream ends of the constriction (see Fig. 3.6. 1-1). 

The pressure and flow are adjusted such that the sample is held at the 
downstream end of the constriction. A suction source is placed 
downstream. A uniform velocity distribution upstream and turbulent 
free flow throughout is provided by antivortex screens in a stilling 
chamber upstream of the constriction. The flow rate is controlled by 
vanes in the upstream volume (Berge et al., 1981) (Rush, 1981). 

An advantage of this levitator is that feedback is not required. 
Disadvantages are that the sample is caused to rotate and that a solid 
sample must be spherical (Rush, 1981). Another disadvantage is that 
the flow may affect the reaction rate (eg. chemical reaction, 
adsorption, etc.) in a gas grain experiment. 

3.6.2 TRIAXIAL SONIC PUMP LEVITATION 

A sonic pump is a loudspeaker driven device that produces a directed 
flow of gas. Six of these may be directed inward in opposing pairs 
along three orthogonal axes and combined with an electro-optical 
detection and feedback system to produce a levitation system. A sonic 
pump is simply a speaker connected to a tube (see Fig. 3. 6. 2-1). 

During one half cycle of the speaker cone, gas enters the tube 
isotropically. During the other half cycle gas leaves the tube with 
directed flow outward. The pulsed nature of the outward flow is 
detectable for only a short distance beyond the end of the tube. The 
levitated object, which is many tube diameters away, experiences a 
steady flow. The advantage of the sonic pump is that it can control 
gas flow with shorter response time than could any mechanical valve 
(Dunn, 1985). 

Sonic pump closures with single orifices are used to levitate large 
spheres, while multiple orifice closures (collimated hole structures) 
are used for smaller objects (Dunn, 1985). 

An advantage of this system is that is it not dependent upon 
temperature (Dunn, 1985). 

3.6.3 TRIAXIAL MOLECULAR BEAM LEVITATION 

For levitation at very low pressures molecular beams may be arranged in 
a triaxial configuration and a vacuum pump may be applied to the 
levitation chamber. Each molecular beam may emanate from a single 
capillary or a collimated hole structure may be used. Because of 
electrostatic charging and polarization forces, it is more difficult to 
levitate conducting particles (Crane et al., 1981) (Rocke, 1981). 
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Figure 3. 6.2-1 Sonic Pump 


3.6.4 PLASMA LEVITATION 


To Study very high temperature (4000 - 10000 K) gas-grain interactions, 
plasma levitation might be used. One or more plasma torches would be 
created by passing the gas through induction coils. The particle would 
be aerodynamical ly levitated in the tail flame of the plasma torch. An 
arrangement for 1 g is shown in Fig. 3. 6. 4-1. (Farnell and Waldie, 
1973). In microgravity, constricted tube or triaxial arrangements 
might be used. 

3.7 ELECTROMAGNETIC LEVITATION 

A conducting particle may be confined at a minimum in an alternating 
magnetic field. A reversed coil configuration (see Fig. 3.7-1 (a)) 
(Wang, 1983) produces a cusp shaped field (see Fig. 3.7-1 (b)) (Frost 
and Change, 1982) which confines the particle at its center. Two 
limitations with this method are that the particle is heated (and 
perhaps melted) and that a minimum particle conductivity of 1000 
(ohm-m)"^ is required (Day and Ray, 1986). 




Figure 3.6.4-1 

Levitation in 1 gin an Induction Plasma Torch 
(Famell and Waldie, 1973) 
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(a) Coil Configuration (Current Direction 
for Peak of Half Cycle) 



Sample 


(b) Field Configuration (Current and Field 
Direction for Peak of Half Cycle) 


Figure 3. 7-1 Simple Configuration for Magnetic Levitation 


4.0 


LEVITATION INDUCED COAGULATION FORCES 


All levitation methods either induce coagulation or otherwise render 
unrealistic the simulation of multiparticle interactions. 

In microgravity levitation forces can be reduced by up to a factor of 
one million. However, the effects of Brownian motion would not be 
reduced at all and for sufficiently small particles the effects of 
Brownian motion must be as important as those of gravity (see Section 
2.3) (aerosols exist in 1 g in the Earth's atmosphere). Obviously, any 
levitation scheme that opposes Brownian motion with sufficient strength 
so as to keep particles away from the walls of the chamber, will 
fundamentally affect the dynamics of small particles in the confined 
cloud and will thereby affect the growth of larger particles. 

4.1 COAGULATION CAUSED BY LIGHT PRESSURE LEVITATION 

Multiple liquid drops have been levitated in a vertical beam in 1 g. 
They become ordered such that the largest is closest to the beam 
focus. Also the upper drops are located in the light intensity maxima 
of the diffraction pattern formed by the lower drops. Up to 20 drops 
have become ordered in a fixed array, which rearranges upon being 
disturbed. Also, irregularly shaped (solid) objects orient themselves 
in the beam (Ashkin and Dziedzic, 1975). 

While the above phenomena may not be strictly defined as coagulation, 
they would probably prevent realistic cloud experiments because they 
would interfere with particle movement and introduce artificial 
anisotropy and diffraction into light scattering. 

Light pressure levitation is the only technique by which "vacusols" 
might be confined. But in the absence of Brownian motion, the ordering 
in fixed arrays would be particularly acute. 
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4.2 


COAGULATION CAUSED BY RADIOMETRIC LEVITATION 


The TEMOl* (doughnut) mode required for radiometric levitation 
(Lewittes and Arnold, 1982) will cause (or at least enhance) 
coagulation on the axis of the beam. Also, highly absorbing carbon 
particles have been observed to assume stable positions in (highly 
localized) energy density minima in a focused TEMOO (gaussian) beam 
(Pluchino, 1983). In either case radiometric levitation will interfere 
with realistic cloud experiments. 

4.3 COAGULATION CAUSED BY ACOUSTIC LEVITATION 

Acoustic levitation will cause acoustic agglomeration, the industrial 
use of which is an intermediate treatment of aerosols that contain 
submicron- and micron-sized particles so that they may be removed by 
conventional cleaning techniques. The high-intensity acoustic field 
causes local velocity fluctuations, which cause the particles to 
collide, adhere, and grow (Tiwary, 1984). This is partly a result of 
the amplitude of acoustic jitter being dependent on particle size 
(Foster, 1985) (see Figs. 4.3-1 thru 4.3-3). 

It is appropriate to estimate the rate of acoustic agglomeration for a 
particle that is acoustically levitated in 10~5 g. Equating Eq. 

2.9-2 to the gravitational force 10~5 mg, using the maximum value of 
the term sin 2kx, and solving for the pressure amplitude yields. 

- ■ 


where Pg is the mass density of the gas and is the mass density of 
the particle. The sound intensity ("a measure of the acoustic energy 
flux") is (Foster, 1985) 


I 



( 4 . 3 - 2 ) 


The sound intensity level in decibels (dB) is (Foster, 1985) 


IL = 101og,o(l/a 


( 4 . 3 - 3 ) 


where Iq = 10~^2 w/m^. For air with a pressure of 765 mb and a 
temperature of 283 K and sound with a frequency of 1000 Hz, Eqs. 4.3-1 
thru 4.3-3 yield a sound intensity level of 125 dB. The rate of 
acoustic agglomeration will be approximated for water droplets with a 
mean radius of 12 um, a radius variance of 1.0, and a total water 
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Figure 4.3-2 Motion of a 10 -fim Droplet (Foster, 1985) 
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Figure 4.3-3 Motion of a 20 -nm Droplet (Foster, 1985) 
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content of 1 g/m^. Letting n(t) be the total number density with the 
initial size distribution, the initial variance is sufficiently small 
that it is approximately true that 


dn 

dt 


-n^ K, 


(4.3-4) 


where K is the collection constant for the initial size distribution. 
[An accurate answer would require the integration of the stochastic 
collection equation in which K is a collection kernel (Foster, 1985).] 
For the stated conditions at 125 dB, K = 0.00006 cm^/s (see Fig. 

4.3-4) (Foster, 1985). The initial density of droplets is 138/cm2. 
Therefore, the initial fractional rate of change of the number density 
of water droplets with the mean initial radius of 12 um is 


1 ^ = .nK = -0.0083/s = -0.50/minute, 
n dt 


Therefore, in one minute, one-half of the droplets initially present 
would have agglomerated. The above calculation is extremely crude. 

One reason for this is that the collection constant used is for 
traveling waves (Foster, 1985) rather than for the standing waves used 
in acoustic levitation. However, standing waves are preferred for 
causing acoustic agglomeration (Barmatz, March 24,1982). 

The kinetic (Brownian) coagulation coefficient (no acoustic 
agglomeration) for particles with a radius of 12 ym is 10x10“^® 
cm’/s (Twomey, 1977), which for an initial density of 138/cm^ 
yields an initial fractional rate of change of -0.36/month. Obviously 
the acoustic agglomeration overwhelms the kinetic coagulation. 

Acoustic agglomeration is not just the same coagulation speeded up. 

One reason is that acoustic agglomeration is highly dependent upon 
particle radius (Foster, 1985). 

Another point is that for organic chemistry experiments that are 
scheduled to last days or weeks or months [The Titan aerosol experiment 
is scheduled to last 90 days (NASA TM-89606, 1987)], acoustic 
agglomeration would have to be taken into account. Questions of the 
following kind should be considered: would the over 75 separate 

compounds comprising the Titan aerosol tholin (Khare, 1984) be formed 
in the same proportion if agglomerated in a few minutes versus three 
months? 
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Figure 4.3-4 

Acoustic Kernels as a Function of Initial Mean Radius (Foster, 1985) 


51 


U.U COAGULATION CAUSED BY ELECTROSTATIC LEVITATION 

For electrostatic levitation with feedback (eg. tetrahedral 
electrostatic positioner) only one particle may be levitated and thus 
coagulation is not defined. 

For levitation with alternating electric fields a circular central 
(alternating voltage) electrode is designed for the levitation of only 
one particle. A rectangular central electrode can accommodate a number 
of particles, which are confined to and free to move along the long 
axis of the rectangle (Straubel and Straubel, 1986). However, since 
these charged particles are mutually repulsive, even this lack of 
coagulation in one dimension does not constitute a realistic simulation. 

4.5 COAGULATION CAUSED BY AERODYNAMIC LEVITATION 

Coagulation is not defined for constricted-tube gas-flow levitation or 
triaxial sonic pump levitation, both of which are designed to levitate 
single particles. 

4.6 COAGULATION CAUSED BY ELECTROMAGNETIC LEVITATION 

Electromagnetic levitation is designed to levitate a single particle. 
With more than one particle present, either the system would fail or 
possibly the particles would coagulate. 
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5.0 


FEASIBILITY OF LONG DURATION EXPERIMENTS 


Single particle long duration experiments are feasible. The 
electrostatic levitation method may require occasional recharging of 
the particle (eg. by a pulse-charging system). The only real limiting 
factor may be the electric power required to maintain the temperature 
of extremely hot or cold chamber walls. 

Since all levitation methods cause coagulation (or introduce some other 
unrealistic effect) it would probably be best to confine cloud 
experiments with chemically inactive walls. The rate of precipitation 
at the walls would then depend on the surface to volume ratio, which is 
minimized for large volumes. Therefore the duration of cloud 
experiments would be limited by the size of the chamber. This suggests 
that two chambers be used: a small one for single particle experiments 
and a very large one for cloud experiments. 

Fractal formation experiments are feasible because the low density of a 
large fractal causes its settling due to residual gravity to be very 
slow. (See Appendix A for an analysis of fractal coagulation.) 
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6.0 


FEASIBILITY OF CLASSES OF EXPERIMENTS 


6.1 


The feasibility will be analyzed of the following classes of 
experiments: 

(a) low velocity collision experiments, 

(b) experiments to simulate high temperature condensation of refractory 
grains, 

(c) cloud coagulation experiments, 

(d) experiments to measure optical properties of clouds, 

(e) Titan aerosol experiments, 

(f) fractal formation experiments, 

(g) coating by coagulation experiments, and 

(h) "vacusol" experiments. 

LOW VELOCITY COLLISION EXPERIMENTS 

The pendelum type collision experiments designed to achieve collision 
velocities down to 10“^ mm/s (based on 10~5 g residual 
acceleration) (Marshall et al., 1985) may not be feasible in general 
because of vibrational accelerations of up to lO"^ g due to astronaut 
motion and equipment vibration. Of course, the motion due to these 
vibrational accelerations could be subtracted out, but then the 
accuracy of the experiment would be limited by that of the 
accelerometers. 


A better procedure might be to place the two particles 1 mm apart in 
the center of the chamber. One would be given a small push toward the 
other and their collision would be recorded with high speed video 
cameras. Vibrations would affect the camera but would have no effect 
on the recorded relative velocity or separation distance of the 
particles. Of course, the residual acceleration would eventually cause 
the particles to hit the wall of the chamber, thus ending the 
experiment. For a chamber of radius 5 cm and a residual acceleration 
of 10“5 g, the particles would hit the wall in 32 s. Thus, for an 
initial separation of 1 mm, the minimum closing velocity of the two 
particles would be 0.03 mra/s. This velocity could be decreased by 
using a larger chamber. However, the size of the chamber is probably 
limited by the field of view of the cameras, although multiple cameras 
might mitigate this problem. 

In the above scheme, vibration may still affect the assignment of 
initial position and velocity to the particles. Therefore a given 
experiment must be performed many times and only these cases having 
initial velocities and positions of interest will be selected for 
analyses. 

6.2 EXPERIMENTS TO SIMULATE HIGH TEMPERATURE CONDENSATION OF REFRACTORY 
GRAINS 

Since all levitation techniques enhance coagulation or cause other 
unrealistic effects, condensation experiments should probably not use 
levitation but be confined by chemically inactive walls. During the 
experiment a sample would be drawn off and an individual grain would be 
electrostatically levitated in a small chamber to study its optical 
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properties. This sample would also provide grains for physical and 
chemical analyses and for the measurement of the particle size 
distribution (see Figure 6.2-1). Once suspended electrostatically the 
grain could continue to grow. The reason for initial nucleation in a 
large chamber is to minimize the surface to volume ratio and thereby 
minimize contamination from the walls. 

Pulsed levitation could be used to position a single grain in the large 
chamber. This could be by laser (radiometric or light pressure) or by 
gas puffs. 

Optical measurements could also be made in the large chamber. Size, 
number densities, and trajectories of individual particles can be 
measured and stored or made available in real-time. One technique is 
the image dissector (Knollenberg , 1979). 

6.3 CLOUD COAGULATION EXPERIMENTS 


These are feasible because they are performed on Earth. The advantage 
of microgravity is that it affords a lower precipitation rate for large 
particles. 

For condensation/coagulation experiments of the atmospheric type, one 
should study the cloud chambers designed for the Atmospheric Cloud 
Physics Laboratory (ACPL), which was designed as a payload for the 
Space Shuttle. A computerized literature search of the Aerospace 
Abstracts yields 42 references to ACPL. 


For general coagulation experiments the appropriate equations must be 
solved to determine if the objectives may be obtained within a 
reasonable time. The time rate of change of total number density is 
easily obtained by solving the following equation (Ludlara, 1980): 

dN 

— = -kN*, (6.3-1) 

at 

where 


k s 8i: r Dps5 X 10“^® cm^/s for air at STP, 


r is the particle radius and Dp is the particle diffusion 
coefficient. The solution is tLudlam, 1980) 


N(t) 


1 

1/No + Kt ' 


(6.3-2) 


where Nq is the initial number density (see Tables 6.3-1 and 6.3-2). 
To obtain the evolution of the detailed distribution, the following 
equation must be solved numerically (Twoney, 1977): 


^ n(x) 
dt 


-n 



K(x,v) n(v) dv 


+ 


V2 


K(x-v,v) n(x-v) n(v) dv, 

0 


(6.3-3) 
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Chamber for Electrostatic 
Levitation and Optical 
Measurements 


Figure 6.2-1 Layout for High-Temperature Condensation of Refractory Grains 






where 

n(x) = number density per unit particle volume for particles 
of volume x (m~°), and 

K(x,v) = coagulation kernel (m^ s~^). 

For coagulation resulting from Brownian motion, K may be written (as a 
function of particle radius) as follows (Twomey, 1977): 

K(rj,r2) = 4 ir (r^ + r2) (D^ + D 2 ). (6.3-4) 

As the Knudsen number approaches zero (Twomey, 1977), 

kT 

D = , and (6.3-5) 

6 n [i, r 

K(r,,r2) = (ti + Ti) • (6.3-6) 

3 \r, r,J 
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Table 6.3-1 


Time t Required to Reduce an Initial Concentration n^ (cm”^) of 
Small Particles of Uniform Size to n^/lO under Ordinary Conditions at 
STP (LUDLAM, 1980) 


(cm ~3) 


t 

lOll 


0.2 sec 

10^ 


5 hr 

105 


2 days 

10^ 


3 wk 

103 


6 mo 


TABLE 6.3-2 


CSiange of Concentration 
One Week, Under Ordinary 

of Small Particles of Uniform Size During 
STP Conditions (LUDLAM, 1980). 

nQ (cm“3) 


n (cm~3) 

10l2 


2.9 X 103 

108 


2.9 X 103 

10^ 


2.2 X 10^ 

2.2 X 103 


1.3 X 103 


58 


Whenever possible, cloud coagulation experiments probably should not 
use levitation but should use confinement by inert walls. In order to 
minimize the effect of wall losses, the chamber volume and the initial 
density should be as large as possible. On Earth smog chambers of 
volume 60 m^ fabricated of FEP Teflon film are used to study aerosol 
formation and even in these large chambers a significant fraction of 
the particles is lost to the walls (McMurry and Grosjean, 1985). 

In addition to particle loss to the walls by Brownian diffusion and 
residual acceleration, turbulent diffusion and electrostatic drift must 
be considered (McMurry and Rader, 1985). The effect of turbulent 
diffusion in microgravity has not been investigated during the present 
study. Electrostatic drift could probably be eliminated by grounding a 
dielectric chamber at numerous points. 

Another consideration is loss to the walls of the gas-phase species 
(according to McMurry and Grosjean, 1985): 

Wall deposition of gas-phase species may also 
affect experimental results. The chemical 
systems that are studied in smog chambers are 
typically complex, involving many gas-phase 
species. The time dependence of any given 
species depends on its rates of formation and 
consumption by other species. Wall-deposition 
rates probably vary from species-to-species and 
thus introduce another removal rate in the 
differential equations that must be solved to 
determine the time history of the species. If 
wall-deposition rates are comparable or large 
compared to rates at which a species is 
consumed by chemical reactions, then wall 
deposition may affect the time history of this 
and other species. 

Since loss rate to the walls of gaseous or particulate species is 
proportional to number density, while rate of particulate coagulation 
is proportional to the square of number density, obviously there can be 
no generalized comparison between these two mechanisms of particle 
removal . 

6.4 EXPERIMENTS TO MEASURE OPTICAL PROPERTIES OF CLOUDS 

Optical properties of individual particles have conveniently been 
measured on Earth via light pressure levitation for transparent 
particles (Lettieri and Preston, 1985) (Ashkin and Dziedzic, 1981) and 
via electrostatic levitation for transparent and nontransparent 
particles (Weiss-Wrana, 1983) (Straubel and Straubel, 1986). These 
experiments could easily be performed in microgravity with an advantage 
being that a liquid droplet would be more nearly spherical when not 
levitated against 1 g. Also, in microgravity the laser levitation 
would cause less heating and evaporation and electrostatic levitation 
would require less charge. 
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The optical measurement techniques listed below were obtained from a 
review paper by Hirleman (Hirleman, 1983). The techniques described 
are all non-intrusive, offering minimal disturbance to the system being 
studied. 

Photography and holography are useful for particles greater than about 
five microns in diameter. They may be automated to determine a mean 
particle diameter, the particle size distribution and concentration, as 
well as specific information such as particle shape (Simmons, 1977) 
(Fleeter et al., 1982) (ICnollenberg, 1977). Using double flash 
photography, in which two closely spaced light pulses are used, the 
particle velocity may be obtained as well (Lennert et al., 1977). 

Pulsed holographic techniques produce three-dimension information on 
particle size and shape (Trolinger, 1980) (Thompson, 1974) (Jones, 

1977). 

A single particle counter is an instrument which shines a laser into a 
relatively small optical sample volume traversed by the particles to be 
measured, and collects scattered light at one or more scattering 
angles. The particle size distribution may be inferred from the 
intensity and/or angular distributions of the scattered light. This 
technique is capable of measuring particle sizes of 0.5 microns and 
above with high specificity (Ungut et al., 1978) (Knollenberg, 1977) 
(Holve and Self, 1979). With a second laser beam, one also has the 
potential for simultaneous velocity measurements (Hirleman, 1978). 

Particle sizing interferometry is a method of particle sizing which 
incorporates laser Doppler velocimetry. In the latter technique, a 
laser beam is split, with the two component beams intersecting a sample 
volume at different angles. A particle passing through the sample 
volume will scatter each beam with a distinct Doppler shift. The 
difference of the scattered frequencies, which contains particle 
velocity information, is obtained by heterodyning the scattered 
signals. In. particle sizing interferometry, the size of the particle 
is then determined, (in an average way, since the particle doesn't 
always traverse the center of the sample volume) from the product of 
the particle velocity, with the time over which the scattered signal 
was received (Farmer, 1972) (Fristium et al., 1973) (Robinson and Chu, 
1975) (C)hu and Robinson, 1977) (Yule et al., 1977) (Farmer, 1978) 
(Bachalo, 1980). 

An instrument for ensemble or multiparticle analysis is the polar 
nephelometer (Hansen and Evans, 1980) (Hansen, 1980), which is capable 
of measuring molecular as well as particulate scattering. In this 
method, light scattering is measured at several angles, and the 
scattering intensities are compared to theoretical values based upon 
Mie theory. In principle, a mean scattering diameter, the particle 
size distribution, and the real and imaginary parts of the index of 
refraction may be evaluated in this manner. A similar technique for 
larger particles, in which the particle size distribution is determined 
from measurements on the forward scattering lobe, is laser diffraction 
particle sizing (Hirleman, 1983) (Chin et al . , 1985) (Swithenbank et 
al., 1977) (Dobbins et al . , 1963) (Roberts and Webb, 1964) (Dieck and 
Roberts, 1970) (Alger, 1979). 
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In diffusion broadening spectroscopy, particle sizes for very small 
particles are determined from their diffusion constants, as measured 
from Doppler shifts due to Brownian motion (Benedek, 1969) (Penner et 
al., 1976) (King et al . , 1982), This technique has the advantage of 
not requiring knowledge of the index of refraction. 

Finally, when the index of refraction and the volume concentration of 
the particles are known, the volume-to-surface-area mean diameter (or 
Sauter Mean Diameter, D 23 ) of the particle distribution may be 
determined by a single transmittance measurement (Dobbins, 1966). 

6.5 TITAN AEROSOL EXPERIMENTS 


As currently proposed the Titan aerosol experiment (McKay, 1986) is not 
feasible. N 2 at the specified pressure of 0.2 mb has a mean free 
path of 3 X 10 “^ m, which is the minimum particle size for acoustic 
levitation (McKay, 1986). For laser light with a wavelength of 6888 
angstroms, a tholin particle of this size is almost opaque (khare, 
1984). For optical levitation against an acceleration of 10"^ g, the 
absorbing power would be 4 x 10~^ W. Assigning an emissivity for 
thermal radiation of 0.1 results in a steady state temperature of 700 
K, which is probably above the starting decomposition temperature for 
tholin. One solution might be to raise the gas pressure, which would 
lower the mean free path and thereby decrease the minimum size for 
acoustic levitation. An obvious problem with this solution is that 
raising the pressure may change the kinetics or even the fundamental 
chemistry of tholin formation. 


Another problem is that when the mean free path of the gas is equal to 
the particle size, the radiometric (photophoretic) force is maximized 
(Ashkin and Dziedzic, 1976). For an absorbing particle such as tholin, 
the radiometric force may exceed the force from light pressure. For an 
incident power P, the photon pressure force is (Lewittes and Arnold, 
1982). 


Fp 



(6.5-1) 


Neglecting the thermal conductivity of the particle and the cooling of 
the particle by radiation, an opaque particle that is much smaller than 
the mean free path experiences a radiometric force of (Lewittes and 
Arnold, 1982). 


F, = 


3v ' 


(6.5-2) 


where v is the molecular velocity. The ratio of forces is 


F, ^ _c_ 
F„ 3v ■ 


(6.5-3) 
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At 300 K the velocity of N 2 molecules is 300 m/s and therefore the 
radiometric force exceeds the force due to photon pressure by a factor 
of 3 X 105. 

A crude calculation shows that the radiometric force can be this large 
even at the low pressure of 0.2 mb. At 0.2 mb and 300 K the number 
density is 4.83 x 10^^ m“3 and the cross section of a 3 x 10“^ m 
diameter particle will encounter 1.02 x 10^^ collisions per second 
with the gas. Assuming a particle mass density of water, the 
levitation force at 10“^ g must be 1.38 x 10“^2 , which for light 

pressure levitation would require a power of 4 x 10"^ W. To remove 
this heat the N 2 molecules would have to leave the lighted surface of 
the particle at a speed of 509 m/s yielding a radiometric force of 
10"^ Nt. Thus, in this very crude calculation, which is based on the 
specified pressure, the radiometric force exceeds the force due to 
photon pressure by a factor of one million. 

The problem with confinement via the radiometric force is that the 
particles would move toward a light minimum, requiring the use of the 
TEMOl* (doughnut) mode. Also, very small particles would be subject to 
negative photophoresis and would move toward the laser, while large 
particles would be subject to positive photophoresis and would move 
away from the laser. Thus as the particle grew the direction of force 
would change. 

A possible way to do this experiment is to let very small tholin 
particles form in a large chamber with no levitation. These would be 
suspended by Brownian motion. A gas sampling tube would then transfer 
a single particle to a small electrostatic levitation chamber (as in 
Figure 6.2-1), where the particle would continue to grow. A possible 
problem is that the artificial charge on the levitated particle may 
interfere with the chemistry of tholin production. 

Another possible method is to use occasional pulses of laser light to 
position a single particle as it continues to grow. 

6.6 FRACTAL FORMATION EXPERIMENTS 

These experiments are feasible and are mentioned in Section 5.0 and 
analyzed in Appendix A. 

6.7 COATING BY COAGULATION EXPERIMENTS 

An experiment has been proposed in which a particle of initial radius 
0.2 ym grows to 0.4 ym by coagulation with 0.1 ym particles and in 
doing so entrains an organic condensate and shields this organic 
material from electromagnetic radiation. This experiment would be 
performed in air at STP and is proposed to last 1-100 hours with a 
constant density of 1-10^/cm^ of 0.1 ym particles. 
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The chemical feasibility of the entrainment of the organic material is 
beyond the scope of the present study. The physical feasibility of the 
coating by coagulation depends on the number density of 0.1 um 
particles. The volumetric rate of growth of the big particle is (from 
Twomey, 1977); 

^ = 4n(D, + D^Hr, + rj - n , (6.7-1) 

dt 3 

where 

(6.7-2) 


(6.7-3) 

n = gas viscosity, and 


D, = 
D, = 


kT 


6 n ri r, ' 
kT 

6 n Ti ' 


n 2 = number density of small particles. 

Since > r 2 , D 2 > Dj and (6.7-1) becomes 

d , 4 „ 3 . 2kT 4 

— (-nr, 3) * r, nj-nr" 


dt 3 
dt 


3 n r 


2 f 


2 kT 

9n 


n2 


(6.7-4) 


Integrating (6.7-4) from initial radius r^j to final radius r^f 
yields 


t 

t 


4 kT rij 
1.16x10’7rij 


(6.7-5) 

(6.7-6) 


Using the maximum number density of lO^/cm^ yields t = 1.16 x 10^ 
s = p days. Therefore the experiment is not feasible for the given 
density, duration, and particle sizes. 

6.8 "VACUSOL" EXPERIMENTS 

"Vacusol" experiments that last more than 30 seconds are not feasible 
because with no gas present, even the smallest particles will be 
accelerated at 10“5 g toward the walls. 
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9.0 


SUMMARY AND RECOMMENDATIONS 


The results of this study may be summarized as follows: 

(a) All levitation techniques either produce artificial coagulation, 
ordering, or other effects that adversely affect cloud 
experiments. Therefore, whenever possible, cloud experiments 
should not use levitation but should be performed in a chamber with 
inactive walls (eg. Teflon). The chamber should be as large as 
possible in order to reduce the surface to volume ratio and thereby 
reduce contamination from the walls. 

(b) Levitation is useful for the study of the optical properties of a 
single particle after it has been nucleated in the large chamber. 

It may be possible to levitate this single particle during 
continued growth. Levitation would be performed in a small 
separate chamber. Electrostatic levitation is a well established 
technique that is probably the most versatile with respect to 
particle size and composition. 

(c) Pendelum-type low velocity collision experiments are not feasible 
because the vibrational accelerations of 10~2 _ io~2 g would 
overwhelm the residual gravity of 10“5 g. However, sufficiently 
low collision velocities could be achieved without a pendelum by 
positioning two particles at the center of the chamber and 1 mm 
apart, giving one a small velocity toward the other, and recording 
their relative motion with a high speed video camera. 

(d) Refractory grain experiments are feasible if small grain sizes are 
acceptable. 

(e) Cloud coagulation experiments are feasible, in general. However, 
the feasibility of any particular experiment depends on having 
sufficiently high initial number density. 

(f) Measurement of optical properties of clouds is feasible. 

(g) As originally proposed, the Tital aerosol experiment is not 
feasible if large particle sizes are required. Proposed 
modifications may enhance its feasibility. 

(h) Fractal growth experiments are feasible. 

(i) The experiment for particle coating by coagulation is not feasible 
for the specified number density, duration, and particle sizes. 

(j) "Vacusol" experiments lasting more than 30 seconds are not feasible. 
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The primary recommendation is that as part of, or prior to any engineering 
feasibility study, the scientific feasibility of each proposed experiment 
should be analyzed in more detail by the appropriate specialists. Depending 
on the experiment, this analysis may require accurate numerical simulations 
(eg. solving the integro-differential coagulation equation) or laboratory work 
on Earth. 
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APPENDIX A 


SELF PRESERVING SOLUTION FOR COAGULATION OF FRACTAL PARTICLES 

The following is a generalization of the derivation (Friedlander and 
Wang, 1966) for ordinary particles (constant total particle volume per 
system volume). The following assumes constant total particle mass per 
system volume. In the following, parts of sentences are borrowed 
without the use of quotations. 

The kinetic equation for coagulation of particles is 


3n (m,t) _ 1_ f” 
at ~ 2 K 


3(m,m-m) n(m,t) n(m-m,t) dm 



3(m,m) n(m,t) n(m,t) dm, 


(A-1) 


where 

m = mass of a particle (kg), 
t = time (s), 

n = number of particles per unit volume in the mass range m to m+dm 
(kg“l m“3), and 

8(&,m) = collision frequency factor between particles of mass ra and m 
(m3/s). 

For coagulation by Brownian motion 

3(m,m) = [r(m) + r(m)] ^ + -i-1 , (A-2) 

3m r(m)J 

where 

k = Boltzmann's constant (l.38xl0“23 j/k), 

T = temperature (K), 

U= fluid viscosity (N-s/m^), and 
r = particle radius (m). 
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For a fractal particle 

m « r, (A-3) 

where z is between 1.7 and 2.5. Using (A-3) in (A-2) yields 


3(m,m) 


2kT 

3|i 


(m’'' + m’^') (m-’'' + m ’'') . 


(A-4) 


The total number of particles per unit volume N(t) and the total mass 
of particles per unit volume 4>> are given by 


N(t) = f 

*’0 


n(m,t) dm, and 


(A-5) 



m n(m,t) dm. 


(A-6) 


The following transformation is introduced: 


n(m,t) 


[N(t)P 

<P 




n = 


N(t)m 


, and 


T = 



(A-7) 


(A-8) 


(A-9) 


I 
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If 6(m,m) is a homogeneous function of its arguments, subs 
(A-7) in (A-1) yields 


{n 


#(n,T) ^{r\.T) dri dff 


^ 2 „ 2T 

an at 


f** ~ 

+ 3(n/n-^ ^(rf,T) $(ri-Ti,T) dff 

-2 ( P(ri,fn dff = 0 . 

« o 


Using (A-7) in (A-5) and (A-6) yields 


f #{ri.T) dri = 1 and 

'>0 

J n #(ri,T) dn = 1 . 


Letting 2 t ^ _ q 

ax 

and using (A-4) in (A-10) yields 


(Friedlander and 


(1 + ab) + (2ab - br|’'‘ - ari ”'') 4){r|) 

dri 

+ 4J(n - 1 + = 0 , 


i tut ion of 


(A-10) 

(A-11) 

(A-1 2) 

Wang, 1966) 

(A-1 3) 
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where 

tj»(n) = replacement for l'(n.T) , which is 

assumed to be independent of T , 


a = I n’'' 4^(n) dn, and 

b = ( n’'' 4»(n) dn- 


(A- 14) 
(A- 15) 


Letting 




(1 + ab) ^-*-^ 
a* 


X(0 


and 


(A- 16) 


n = 


a" 

(1 + ab)‘ 




(A-13) becomes 


? ^ + [2a-a(1 -a)C - Cl X(^;) 

d^ 




(A-17) 


(A-18) 


where 


ab 


a 


(1 + ab) ’ 


(A- 19) 



In terms of x(C), Eqs. (A-ll), (A-12), (A-14), and (A-15) become 


[ x(^:) d^; = 1 - a , 


i' = ?77bF- 


fVx(2;)d^ = 1,and 


( d^ = an - a) 


(A-20) 


(A-21) 


(A-22) 


(A-23) 


For the lower end of the spectrum, the integral term in (A-18) may be 
neglected: 

Z, + I2a - an - a) C - X(0 = 0 , (A-24) 

dq 


a solution to which is 


X(0 = exp Iz an - a) C - z C’''l . 


(A-25) 


where C^, is an integration constant. For the upper end of the 
spectrum, the second term in (A-18) is negligible: 



(A-26) 
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N| CM 


A solution to which is 


x(^) 


— — 

r+T7TT^v^7TT~vI) ' 


(A-27) 


Where C 2 » is an integration constant and B is the beta function. 
Using the substitutions 


^ = e*, and 
y(x) = x(0. 

Eq. (A-18) may be written as 

[2a - a(1 - a)e"'' - e *'*] y(x) 

[ y{^ y[ln (e* - e*)] 

J -OO 

1 + ’'‘j e*" dx = 0 . (A-30) 

Using the substitution 

t = - (X - x) , 

the integral in (A-30) may be written as 


dy(x) ^ 
dx 

+ 


(A-28) 

(A-29) 


6* i y(x-^t) y[x + ln(1 -e’^‘) [e’^ + (1 - e’^)’''] 

•'o 2 

IL^ t 

. e 2 dt. 


(A-32) 
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Eq. (A-30) was solved for z=2 (m ® r^) over the range in which y(x) 
is significant: -4.4 £x £6.0.. The starting value y(-4.4) was given 
by (A-25) with and a guessed until Eqs. (A-20), (A-22), and 

(A-23) were satisfied. Then C 2 was chosen to match the upper end 
solution. Using crude stepwise integration (see Table A-l) on a very 
slow desk top computer yielded the following very approximate values: 

Cl = 0.4 

a = 0-6 

a = 0.90 

b = 0.42 

C 2 = 0.2 


The fit for the upper end solution was very poor and the correct value 
of C 2 may be an order of magnitude higher or lower than 0.2. 

However, this value of C 2 is sufficiently accurate to give order of 
magnitude times for formation of large particles. 

The time to form a unit density of particles larger than mass mi, 
will now be computed. Using the above values of a, b, and C 2 in 
(A-16), (A-17), and (A-27) yields 


4^(n) = 0.25 


(A-33) 


Using (A-7), (A-8), and (A-33) yields 


n(m,t) 


[N(t)P 

<P 


0.25 exp 


N(t)m ' 


(A-34) 


At time t, the density of particles with mass greater than mi, is 


N,(t) 



n(m,t) dm = 0.54 N(t) exp 


-0.47 


N(t) m, ' 
< 1 > 


(A-35) 
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Table A- 1 Program for Solution of Equations A-30 Thru A-3 2 

1 ! CoagS 

5 DIM YU000) 

10 Z=2 
20 filfa=.60 
30 Cl=.4 
40 Xtnax=5.4 
50 Xmin=-4.4 
60 Nx=400 

70 Dx=<Xmax-Xmi n)/Nx 
75 Dt®2/2»Dx 
80 2etal=EXP<Xtnin> 

90 Expl=EXP<2*fll fa)*2etal^< l/ZJ-Z^Zetal-'C-l/Z) ) 

100 Y<l)=Cl/Zetal^<2*Fi1fa)»Expl 

105 1=0 

110 X=Xmin-Dx 

120 X=X+Dx 

130 1=1+1 

140 Terml=<2+nifa-fl1fa*<l-fllfa)*EXP<X/2)-EXP<-X/Z))#Y<I) 
145 Sum=0 

147 IF X=Xmin THEN GOTO 250 

150 Ttnax=<X-Xmi n)*2/2 

160 Tmin=-<2/2)*L0G<l-EXP<Xmin-X)) 

165 IF Tmin>Tmax-Dt THEN GOTO 250 
170 T=Tmin-Dt 
190 T=T+Dt 

200 Il=INT<<X-Z/2+T-Xmin)/Dx+l. 1) 

210 I2=INT<<X+L0G<l-EXP<-<Z/'2)»T))-Xmin)/Dx+l. 1) 

220 Term2=EXP<-T/2) + <l-EXP<-2*T/'2))''<l/Z) 

230 Sum=Sum+Y< 1 1 >*Y< I2)*Term2*EXP< < l-Z)*T/2) 

240 IF T<Tmax-l. l»Dt THEN GOTO 190 

250 Y<I + l)=Y<I)-Dx*<Tertnl+Dt*<Z/'2)*EXP<X)*Sutn) 

260 PRINT •• X=";X + Dx,“ Y=“;Y<I + 1> 

270 IF X<Xmax-l. l*Dx THEN GOTO 120 

280 Suml=0 

290 Sum2=0 

300 Sum3=0 

310 Sum4=0 

315 X=Xmin-Dx 

320 FOR 1=1 TO Nx 

330 X=X+Dx 

340 Zeta»EXP<X) 

350 Dzeta=Zeta«Dx 

360 Suml»Suml+Y< I )*D 2 eta 

370 Sum2=Sum2+Zeta+Y< I >»Dzeta 

380 Sum3=Sum3+Zeta''< 1/'Z>*Y< I )*Dzeta 

390 Sum4=Sum4+Zeta'^<-l/'Z)*Y< I ) + Dzeta 

400 NEXT I 

410 PRINT " Input Cl = ";C1 

420 PRINT •• Input Alfa = ";fnfa 

430 Kl=Suml 

440 K2=Sum2 

450 K3=Sum3 

460 K4=Sum4 

470 PRINT '■ K1 = ”;K1 

480 PRINT “ K2 = " ; K2 

490 PRINT K3 = ";K3 

500 PRINT " K4 = ";K4 

510 Ck 1=1-A1 fa 

520 Ck4=fll fa+< 1-ftl fa> 

530 PRINT " 1 - Rifa = ” ; Ck 1 
540 PRINT " fll fa*< 1-Rl fa) = ";Ck4 
610 STOP 
620 END 




The total density of particles is (Ludlam, 1980) 


N(t) 


1 

1/N, + Kt ' 


(A-36) 


where 

N, = N{0), and 

K » SxlO ’® cmVs = 5x10-’® m=>/s. 


Since Nq > 10®/cm^ = 10^^/m^ and t 
to be considered), 

N(t) « (Kt)-’. 

The mass m^, is given by 

m, = m. 0 = , 

while the total mass density is 

<p = N„m,. 


> 1 hour (for the cases 

(A-37) 

(A-38) 

(A-39) 


Using (A-37) through (A-39) in (A-35) yields 



(A-40) 


Setting Nj(t) * 1/V where V = 0.01 m^ = (8.5 inches)^ = useable 

volume and using r^ = 30 nm = 3 xl 0~8 m and r^ = 1 mm = lO”^ m 

yields the second column of Table A-2. The particle volume fraction is 

fvoi (t) = f V n(m,t) dm = 9.2 (A-41) 

•*0 

and this is listed in the third column of Table A-2 corresponding to 
the times of the second column. It is apparent from Table A-2 that for 
initial number density less than lO^^Vcm^ the time is unreasonably 
long, while for the initial number density greater than lO^^/cm^ 
the final volume of particles exceeds the system volume. For an 
initial density of 2.75 g/cm^, for a particle of initial radius 30 
nm, a particle of final radius 1 mm will have a final density of 
8.25x10-5 g/cm^. Using this density and a chamber radius of 
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Table A-2 Time to Create One Particle Per 0.01 ra^ with Radius Greater 
than 1 mm. 

Initial Density (cm~-) Time Final Volume Fraction 


lo8 

35 years 

1.8x10-^ 

109 

2.7 years 

1.6x10-3 

lOlO 

81 days 

1.5x10-2 

10^1 

7 days 

0.14 

10^2 

14 hours 

1.2 

10l3 

1.3 hours 

12 

lOl^ 

7 minutes 

114 


Conditions; 

initial radius = 30 nm 
gas: standard atmosphere 


82 


0.5 m, Eq. 2.3-3 yields a time of 57 days for a 1 mm radius particle to 
reach the wall due to a gravitational acceleration of 10“^ g. 

Therefore, residual gravity is not a problem. Also, from Table 2.3-1 
it is apparent that Brownian motion is not a problem even for the 
initial particles. 
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APPENDIX B 


COMPARISON OF ELECTROSTATIC LEVITATION OF A PARTICLE HAVING A NET 
CHARGE WITH THAT OF A DIELECTRIC PARTICLE HAVING NO NET CHARGE 


A dielectric particle with no net charge will be levitated in the field 
of a spherical conductor of radius An outer sphere (not 

necessarily conducting) of radius rQ^^er will serve as the outer wall 
of the chamber. The center of the dielectric particle will be at 

•^mid’ 


^mid 


^Jnner 


(B-1) 


The particle is a cube with side Ar , positioned such that the 

normal to one side is in the radial direction. An electric field 
normal to a thin planar dieletric induces a surface charge per unit 
area qi/A on one face and a surface charge per unit area -qi/A on 
the other face such that (Halliday and Resnick, 1962) 

D = S, E + ^, (B-2) 

A 


where 

D = electric displacement based only on free charges (neglecting 
qj) (C/m^), 

= permittivity of free space (8.854x10"^^ F/m), and 
E = electric field strength (V/m). 

For an isotropic media (Halliday and Resnick, 1962) 

D = k!e, E, (B-3) 


where < is the dielectric constant and is equal to unity for a vacuum. 
Using (B-3) in (B-2) yields 


9l 

A 


D 



(B-4) 
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Eq. (8-4) has been derived for an infinite planar slab of dielectric. 
The crude approximation is now made that Eq. (B-4) applies to the cubic 
particle : 


(Ar)^ V k/ 

q, = (Ar)- D.. ■ 


where is D evaluated at Since for a positive potential 

on the inner sphere the inner surface of the particle has a charge 
-qj^ < o and the outer surface has a charge qi > o i the net force 
on the particle is 


Fe = -q, - Ar/2) + q, E(r„« + Ar/2). 


(B-6) 


It will be assumed that each charge resides just inside the surface of 
the particle such that 


E(r^ - Ar/2) 


D(rn^ - Ar/2) ^ Q 

K e, 4n K e, {r„^ - Ar/2)^ ' 


(B-7) 


E(r^ + Ar/2) 


D(r^ + Ar/2) _ q 

K e„ 4n K Eo (r„i, + Ar/2)2 ' 


(B-8) 


where 1 and q is the charge on the spherical conductor. 
Also 


^ , and (B-9) 

4n rj^ 


q = 4ne, (B-10) 
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where ^inner voltage on the spherical conductor. Using (B-5), 

(B-7), (B-8), (B-9), and (B-10) in (B-6) yields 


F. = 


So (Ar)^ 

r„J 


1 


1 


- Ar/2)^ + Ar/2)^ 


(B-11) 


The term in brackets in (B-11) is approximately equal to 2 Ar/r„y^ 
such that 


F. = - 


2e, (K- 1) (Ar)^ 


r, 


mid 


(B-12) 


It will be assumed that an acceleration of magnitude 10~^g is 
directed radially inward; 


F, = -{Ar)^p 10»g, 


(B-13) 


where p is the density of the particle. Equating (B-12) and (B-13) 
yields 


^ p 1 0-»g ^ 

V 2(K - 1 ) e. 


(B-14) 


A typical dielectric constant is that of quartz (E perpendicular to 
optic axis): 


K = 4.34. 
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Using a density of 


p = 1 g/cm® (not that of quartz) 


and radii of 

*^inner ” 

^outer ~ cm, and 
^mid “ 

yields Vj^ngj. = 1.54x10^ volts. It may be noted that the size of 
the particle does not enter into Eq. (B-14). 

To levitate a spherical particle charged with n^ electrons in the 
arrangement described above requires that 





10*g 


t|now ® 


t 


(B-15) 


where Rp is the radius of the particle. Eq. (B-15) is evaluated in 
Table B-1. 
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Table B-1 Conditions for Levitation of a Charged Particle 


Particle Radius 
(m) 

1 

10 

100 

1000 


Number of electrons times voltage (V) 
on inner spherical wall 
0.577 
577 

5.77x105 

5.77x108 


Conditions: 

Radius of inner wall = 10 cm 

Radial position of particle = 15 cm from center 
Density of particle = 1 g/cm^ 

Residual gravity = 10“^ g 



Levitation at lO'^g via the dielectric effect is probably not 
practical because it requires a voltage on the inner wall of 1.54 x 
10^ volts. On the other hand, levitating a particle of radius 1 um 
charged with one electron (With the electron at the center of the 
particle, the surface voltage would be 1.44x10“^ V.) requires a 
voltage on the inner wall of only 0.577 volts. Also, levitating a 
particle of radius 1 mm charged to a surface voltage of 10 volts, 
requires a voltage on the inner wall of only 83.3 volts. 
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APPENDIX C 


AMPLITUDE OF MOTION DUE TO BROWNIAN MOTION 
WITHIN AN ACOUSTIC POTENTIAL WELL 


For a one dimensional standing wave the acoustic radiation pressure 
force is given by (see Eq. 2.9-2) 

F = K sin 2kx, (C-1) 

where 

K = — k (C-2) 

6 p 


Since only an order of magnitude answer is desired, Eq. (C-l) will be 
applied to a spherical potential well, with x being the radial 
coordinate and x=o being the center of the well. The translational 
energy of a particle due to Brownian motion in a gas is 



where IC 3 is Boltzmann's constant (l.38xl0~2^ J/K). Therefore the 
amplitude of the motion Ax satisfies 


J Ax 

F dx'. 

o 


(C-4) 


Using (C-l) in (C-4), integrating and then using the approximation 


COS z 



(C-5) 


yields 


Ax 


3k, T 

2 K k 


(C-6) 
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As a check, an approximation to this result will be derived by equating 
the time for diffusing from the center of the well to Ax to the 
time for the force F to return the particle from Ax to the center of 
the well. From Eq. (2.3-1) the time to diffuse to Ax is 

t . (C-7) 

2k, TB 


The time for the force F to return the particle from Ax to the center 
of the well is 


t 



(C-8) 


where 


V = FB, 


(C-9) 


F « K sin 2k -^ = K sin kAx » K kAx. (C-10) 

Using (C-7) thru (C-10) yields 


Ax 



(C-11) 


which is greater than (C-6) by a factor of only /3/2 = /4/3 

= 1.15, which implies an error of only 15X. Balancing the maximum 
value of the force (when 2kx = ir/2) against an acceleration of 10“5g 
yields 


K 



pp 1 0 ®g. 


(C-12) 


where p is the mass density of the particle. The wavenumber is given 
by 

. _ 2n _ 2n_f (C-13) 

X c 


C 
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where 


f = frequency (Hz), 
c = speed of sound (m/s). 



(C-14) 


Y - ratio of specific heats (1.40 for diatomic gases), and 
m = mass of gas molecule (kg). 

Using (C-12) and (C-13) in (C-6) yields 


Ax 


9 kaT c 


16n^ pp 10 ®g f 


(C-15) 


Therefore the largest Ax is produced by the smallest particle radius 
a. M. Barmatz of JPL (private communication via G. Fogleraan of RCA, 
July, 1987) has estimated that streaming limits the minimum particle 
radius to 1% of the wavelength: 

a = 0.01 K and 

a = 0.01 c/f. (C-16) 


Using (C-14) thru (C-16) yields 


Ax 


750 / m ^ 

n ^ypplO-^g 


(C-17) 


For nitrogen gas: 


7 = 1.40, and 
m = 4.65x10^«kg. 
For 


pp = 1 g/cm\ 

Ax = 1.39x10 ’° f. (C-18) 
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Equations (C-16) and (C-18) are evaluated in Table C-1, from which it 
is evident that for practical frequencies (« MHz) the amplitude is 
less than the radius of the particle. Therefore acoustic levitation in 
residual gravity prevents Brownian motion. 


93 


Table C-1 Amplitude of Brownian Motian in an Acoustic Potential Well 


f (Hz) 

a(m) 

Ax (m) 

Ax/a 

10^ 

3.53x10-3 

1.39x10-'^ 

3.94x10-5 

10^ 

3.53x10-*^ 

1.39x10-6 

3.94x10-3 

105 

3.53x10-5 

1.39x10-5 

3.94x10-1 

iq6 

3.53x10-6 

1.39x10-'^ 

3.94xl0l 
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